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"-^ I Abstract. This talk is a brief review of the different methods of galaxy 

fr ■ cluster mass estimation. The determination of galaxy cluster mass is of 

great importance since it is directly linked to the well- known problem 
of dark matter in the Universe and to the cluster baryon content. X-ray 
observations from satellites have enabled a better understanding of the 
^ ■ physics occuring inside clusters, their matter content as well as a detailed 

^D , description of their structure. In addition, the discovery of giant gravita- 

Jil^ ■ tional arcs and the lensing properties of clusters of galaxies represent the 

pg . most exciting events in cosmology and have led to many new results on 

(^ ' mass distribution. In my talk, I will review some recent results concerning 

C^ . the mass determination in clusters of galaxies. 

Q ■ 1. Introduction 

c/3 . Clusters of galaxies are the most extended gravitationally bound systems. They 

^ ' provide an ideal tool for cosmologists to study the formation and evolution of 

the structures of the Universe. They present strong evidence for the presence of 
large amounts of dark matter. Therefore it is essential to determine in a very 
accurate way their gravitational masses to better constrain the still unknown 
a I cosmological density parameter Oq. 

Historically, the evidence of the presence of a huge missing mass was derived from 
the application of the standard virial theorem (Zwicky 1933), which is based on 
the assumption that mass follows the light distribution, but this assumption 
has not yet been confirmed. In this talk I will show that the total cluster mass 
depends on the relative distribution of the visible and invisible components and 
I will discuss the accuracy of the masses derived under the mass-follows-light 
assumption. 

Clusters of galaxies, are also strong X-ray emitters. Since the discovery of the hot 
diffuse gas responsible for X-ray emission, astronomers have started to use X-ray 
observations to constrain cluster masses. Methods based on such observations 
have several advantages compared to optical methods. However, it is not yet 
clear how accurate the standard methods such as the hydrostatic /3-model are. 
ROSAT observations of the Coma cluster have led to a large fraction of baryons 
in contradiction with the standard Big Bang nucleosynthesis predictions. This 
baryon catastrophy has several implications for cosmology in particular on the 
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value of the density of the Universe, Qq . 

Finally, the detection of gravitational lensing in clusters of galaxies has provided 
astronomers with the most powerfull tool for mapping the mass distribution. 
The mass estimates using the lensing method are in general in good agreement 
with the optically derived masses while the X-ray method has systematically 
underestimated cluster masses by a factor 2-3. I will first describe briefly the 
observational properties of clusters of galaxies, then I will review different meth- 
ods which are usually used to estimate their masses and discuss their reliability. 
In this paper, I will adopt the value of Hq = 50 km s~^ Mpc 



2. Observational properties 

2.1. Optically 

Optically clusters of galaxies appear as large concentration of galaxies in a small 
volume. A typical cluster has several hundred ~ 1000 of galaxies, which are 
mainly ellipticals and SOs in irregular clusters. The typical scale radius is about 
IMpc. The distribution of these galaxies has most traditionally been fit by an 
isothermal gravitational sphere which has the approximate analytical form given 
by King's model 

ngaiir)(x[l+r/R,]-\ (1) 

where r is the projected radius and Re ~ O.SMpc is the core radius in a typical 
cluster. 

The radial velocities of the cluster members in a well-relaxed cluster are dis- 
tributed according to a Gaussian distribution. 

N{vr)dvrC<: ex.p{—Vr/2a'^), (2) 

where a^ =< {vr— < w >)^ > is the line- of- sight velocity dispersion 
Merritt (1994), has shown that the mass distribution can be constrained from 
an analysis of the shape of radial velocity histograms, but his method requires 
a large number of measured radial velocities. Indeed, a redshift survey of rich 
clusters of galaxies has typically ~ 50 velocity measurements per cluster, making 
this method unusable except in the case of the well-studied rich cluster Coma 
which has ~ 600 measured radial velocities, but still, there is the problem of 
substructures. 

2.2. X-ray emission 

The X-ray emission from clusters of galaxies is mainly due to hot and diffuse 
intra-cluster gas with T^ ~ 10^ — 10* K and a central density of nx{0) ~ 
10~^cm~'^ (see the excellent review by Sarazin 1986) 

This hot intracluster gas is the main baryon component of clusters of galaxies: 
its mass is several times that of the stellar mass Mg ~ 5 — 7M^, (David et al. 
1994). It represents a large fraction of the total mass (visible-|-dark matter) and 
can reach values of 30% of the total binding mass (Bohringer, 1994). This gas 
radiates by thermal bremsstrahlung emission 
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For very hot gas the spectrum is dominated by the continuum and the only hne 
which is detected in this continuum is the iron hne. At cooler temperatures 
however, some heavy element emission lines such as O, Si, S, Ar and Ca start 
to appear. 

The detection in the X-ray spectra of the iron K-line at 6 kev has shown that 
the gas has been enriched in metals. These metals have been processed into 
cluster galaxies and ejected into the ICM through SN driven winds or outflows, 
providing evidence of a non-primordial origin of part of the gas. The typical 
abundances are about 1/3 -1/2 solar (Mushotzky 1996). What is the quantity 
of the ejected gas ? And what type of galaxies enriched the ICM? All these 
questions are still open (Arnaud 1994). For the mechanism of metal enrich- 
ment of the ICM, it is now well accepted that supernovae are reponsible for 
the injection into the ICM of the heavy elements processed into stars but we do 
not yet understand the relative importance of both types (Matteuci this school) . 



2.3. The Baryon Catastrophy 

Standard Big Bang nucleosynthesis predictions of the primordial abundances 
place tight limits on the present day baryon density in the Universe, 

0.04 < Ofe^^^/i^ < 0.05 

Walker et al. (1991). This is only a small fraction of the critical closure density 
of the Universe. 

White et al. (1993) have noted that hot gas in the Coma cluster contributes 
~ 15% of the total mass within the Abell radius. Thanks to the wide field of 
view and high sensitivity of the ROSAT satellite, it has been possible to reliably 
measure the baryon fraction of the Coma cluster to an even much larger radius 
~ 4Mpc (Briel, Henry & Bohringer, 1992), where this fraction reaches the value 
of 30%. If dark matter is distributed similarly to the X-ray gas, the conser- 
vative value of the gas fraction in Coma cluster /{, ~ 15% leads to ilf, ~ 0.15, 
which is ~ 3 times the universal ^i,^^^ value. Previous X-ray analyses of galaxy 
clusters with the Einstein and EXOSAT observatories have already found high 
baryon fractions, but the authors have not emphasized the implications of such 
quantities of baryons. More recently, compilations of X-ray cluster data and 
their analysis by White & Fabian (1995, hereafter WF) and David, Jones & 
Forman (1995, hereafter DJF), have led to the same conclusion, showing that 
the problem of baryon overdensity is common in clusters of galaxies. What are 
the cosmological implications of this result? The most obvious one is that Oq 
is less than one. Indeed, one way to reconcile the baryon fraction from clus- 
ter analyses (/{, ~ 15%) with the primordial nucleosynthesis prediction, is that 
0.26 < ^oh^ < 0.33. That means that the Universe is open. Recent measure- 
ments of the primordial deuterium abundance D/H from quasar absorption line 
spectra have produced two different values, a low value fi^/i^QQ = 6.2 it 0.810^'^ 
(Rugers t Hogan 1996) and a high value ^bhjoo = 0.024±0.006 (Tytler, Fan 
& Buries 1996). If one accepts the higher D/H value and accounts for baryons 



within clusters of galaxies then Qgh^''^'-^ 0.6. The J7o = 1 universe can be res- 
cued if one believes either in a low Hubble constant (Hq < 40 Figure 1, Bartlett 
et al. 1995, Lineweaver this volume), or in a non-zero cosmological constant 
such that Q.O = ^matter + ^Ai where Q.a = A/3H'^ but still, this is not consistent 
with dynamical evidence of large Oq. 
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Figure 1. This plot is from Lineweaver et al. (1997). It gives the region of 
the h- fib plane allowed by various constrains:BBN, cluster baryonic fraction 
(White et al), and the shape parameter T from galaxy and cluster scale 
density Ructuations. A low h value is preferred 

Other possible solutions are: 

1.- The calculations of standard primordial nucleosynthesis are incorrect. 

2.- The X-ray gas is more concentrated with respect to the dark matter but 
White et aZ.(1993) have shown that gravitational and dissipative effects during 
cluster formation cannot account for such baryon overdensity. 

3.- The intracluster gas is multi-phase, but no model has been proposed to 
explain such a clumpy configuration of the gas. 

4.- There is a problem with the mass estimates. This solution will be dis- 
cussed in next sections. 



Finally, one may ask another interesting question. Is the baryon fraction the 
same in all clusters ? In the standard picture of cluster formation driven solely by 
gravitational instability and where cluster evolution is entirely self-similar, the 
expected baryon fraction should be constant, because no segregation between 
the gas and dark matter has occured. However, if we gather all the derived 
baryon fractions in the literature and compare them, the answer is clearly NO. 
For example the derived mean values of WF and DJF samples are different, 
15% for the former and 20% for the later (but see Evrard 1997). More recently, 
Lowenstein &: Mushotzky (1996) have shown evidence of variations in baryon 
fraction from their analysis of two poor Abell clusters, A1060 and AWM 7, 
using the most recent X-ray observations from ROSAT and ASCA (Figure 2). 



Such variation in baryon fraction from cluster to cluster requires some process 
in addition to gravity, like feedback mechanisms or some other non-gravitational 
effects as suggested by DJF, but there are no theoretical arguments justifying 
such ideas (White et al. 1993). 
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Figure 2. This plot is taken from Lowcnstcin & Mushotzky (1996). It 
gives the enclosed baryon fraction versus radius in units of the virial radius 
in A1060 and AWM7. The dotted (solid) line represents the best-fit mass 
model, the dot-dashed (dashed) lines the most compact and diffuse models 
for A1060 (AWM 7). The compact models lie below the best fits. 



3. The dark matter problem 

Galaxies and hot gas are only a small part of the total cluster mass. The dom- 
inant component is the dark matter. Zwicky in 1933 then Smith (1936) have 
shown that the virial mass exceeds by a large factor the luminous mass. This 
led them to invoke for the first time, the problem of missing mass. To quantify 
the amount of dark matter, we usually calculate the mass-to-light ratio (MLR). 
The mean value found for rich clusters using blue luminosties is M/Lb ~ 3001i. 
In this unit the M/Lb required to close the Universe is ~ 1200ii. Therefore, if 
we assume that clusters of galaxies are good tracers of the whole Universe, then 
ilo = 0.2 — 0.3. So, if one believes in a matter dominated flat cosmological model 
ilo = 1) then where is the missing mass? Still now we do not know how the dark 
matter is distributed relatively to the visible matter. Cluster mass determina- 
tions using optical observations are based on the assumption that mass follows 
the light distribution which is just an assumption not yet confirmed. What are 
the predictions on the mass distribution? Cosmological theories predict that 
dark matter is more diffusely distributed than galaxies. West and Richstone us- 
ing N-body simulations have indeed confirmed this behaviour (West &: Richstone 



1988). Furtheremore Hughes (1989), using X-ray observations of the Coma clus- 
ter, has shown that models where dark matter parallels the distribution of hot 
gas are ruled out by the data. With the improvement of weak lensing analysis, 
one may hope that this question will be answered more precisely, in the near 
future. 



4. Mass determinations 

In this section, I will review the methods used to estimate clusters masses and 
will discuss their validity. 



4.1. Optical methods 

The Virial Theorem method : 

Early estimates of cluster masses (before X-ray observations became possible) 

were based on the application of the Virial theorem. If one assumes that clusters 

of galaxies are bound and self-gravitating systems then the virial mass is given 

by: 

M.~3^, (4) 

where o"^ and Rg are evaluated from the radial velocity distribution (2) and 
the projected spatial distribution of a fair sample of galaxies. As we have seen 
the naive application of equation(4) leads to large amounts of dark matter. 
Therefore the question we want to address here is: How secure are the virial 
mass estimates? Projection effects, contamination by foreground galaxies and 
anisotropy of the velocity distribution may introduce uncertainties into the de- 
termination of the mass. But they are small effects and can not explain such 
large virial masses. Several observations at both optical and X-rays wavelengths 
provide convincing evidence of the presence of substructure in a large sample of 
clusters (Baier 1983, Bird, 1994, Mohr et al. 1993). X-ray imaging observations 
with the Einstein satellite first revealed such complex structure (Forman et al. 
1981) in contrast to the smooth shape assumed in previous studies. Even clus- 
ters that exhibit a fairly smooth and apparently well-relaxed configuration , like 
the Coma cluster, have been found to contain substructure (Fitchett & Webster 
1987, Mellier et al. 1988) with a large subcluster centered on NGC 4839 that 
appears to be falling into the Coma cluster. If this subclustering is not correctly 
taken into account, this would introduce large uncertainties in the dynamical 
mass. A substructure with 10% of the mass can introduce an underestimation 
of 40% on the MLR. 

However, the most serious problem of using the virial theorem comes from the 
fact that we do not know how the dark matter is distributed. Indeed, the 
application of the standard virial theorem assumes that mass follows the light 
distribution. What happens when this assumption is relaxed? It has been shown 
(Sadat 1995) that in this case the standard application of the virial theorem in- 
troduces a bias on the cluster masses, and this bias (/U = {M/L)(jyn/{M/L)true) 
depends strongly and in a non-linear way on the relative concentrations of the 



visible and invisible components. It is found that the cluster mass is over-(under) 
estimated if the dark matter is more (less) concentrated by an amount 
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[1 + 2CcRtruc + C'A-^true] 
[1 + CyRtruc] 



(5) 



where -Rtruc is the true ratio of the masses MDu/Mgi^i and C^,, C\, Cc are the 
relative concentrations of the 2 components. As an illustration of this effect 
we have plotted in Figure 3 the bias /j, versus i?truc in the case where the dark 
matter is less concentrated than the galaxies. 




Figure 3. The bias ^ versus the true ratio R = Rtme of dark/luminous 
mass in the case where the mass is more diffusely distibutcd than galaxies. 
The three curves correpond to different concentrations 

Figure 3 shows that for a ratio of, say, i?truc~ 30, the virial theorem leads to a 
dynamical ratio ~ 3 to 7 times lower. In this case the virial mass determination 
underestimates the true mass, while in the case where the dark matter is more 
concentrated (Figure 4) the dynamical ratio reaches the value of ^ ~ 100 which 
is 3 times higher than the true value ~ 30!. The true mass is overestimated. 
Note, that in the mass- follows- light case, Cc = Cx = C^ = 1 and fi= 1 + 
Rtruc- If one defines a new quantity Rdyn = A* ~ 1 which measures the virial 
estimated ratio of dark matter mass to visible mass, one can see that in this 
case Rdyn = Rtruc and therefore that the virial mass is equal to the true mass. 
It seems clear from this analysis, that as long as we do not know anything on 
the distribution of dark versus visible matter, one has to be "sceptical" about 
the masses derived from the "virial" method. 



Kinematic method If the system is in equilibrium, one can use the equation of 
stellar hydrodynamics to derive the mass 
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Figure 4. Logarithm of the true ratio R versus the bias /x in the case where 
the mass is more concentrated than galaxies. 



For isotropic orbits, 0"^ =crt, there would be a unique solution. Unfortunately 
the orbits of the galaxies are poorly known. Therefore, we have to solve this 
equation with three unknown quantities: crr{r), crt{r) and M(r). Generally, it 
is assumed that either /3(r) = 1 — at{r)'^/crr{r)'^ and M(r) are known functions 
and then derive Ugai and (Tr(r) from eq. (6) which are consistent with obser- 
vations. Unfortunately, the observed velocity dispersion profiles of clusters of 
galaxies are poorly known and can not put strong constraints on the mass. In- 
deed, even for the best studied Coma cluster, Merritt (1994) has shown that the 
observed velocity dispersion profile of this cluster is consistent with several mass 
distributions. 



4.2. The hydrostatic isothermal /3-model 

Problems encountered with optical methods like the shapes of galaxy orbits, the 
small number of galaxies in a cluster, effects of contamination and projection 
can be avoided by using the observations of the hot X-ray emitting gas. The 
gas can be treated as an isotropic fluid, since the elastic collision times for ions 
and electrons are much shorter than the timescales for cooling and heating. The 
timescale required for a sound wave in the intracluster gas to cross a cluster is 
given by tx = l.3[Tg,s/W^r'^^[R/lMpc] Gyr. 

Furtheremore, since this time is shorter than the dynamical time of the cluster 
(~ 10 Gyr), the gas can be assumed to be in hydrostatic equilibrium with the 
cluster potential (Sarazin 1986). 

Under the assumption of spherical symmetry the equation of hydrostatic equilib- 
rium (balance between the pressure and the gravitational forces) can be solved 
for the mass interior to r, M(r): 
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where Tgas{T:) and pgasi^) are the temperature and the gas density profiles, k is 
Boltzmann's constant, and fi rrip is the mean molecular weight of the gas. In 
principle, the knowledge of Tgas and pgas from the observations, directly yields 
the actual mass distribution M(r). This method has several advantages over the 
optical approach. The gas is isotropic, there are no contamination effects and 
the most important advantage is that the mass distribution is derived directly 
without any assumption about the dark matter distribution as is the case with 
the optical method. 

The sad point is that one must recover three dimensional profiles from pro- 
jected profiles. For the temperature information, this requires the measurement 
of Tgas{T:) which is still very difficult to obtain, even with the ASCA satellite. In 
practice, we assume that the gas is isothermal at a mean temperature Tx- Nu- 
merical simulations (Evrard 1996) and recent ASCA results (Ikebe et al. 1994) 
seem to support this assumption at least out to a radius of 1.5 Mpc. If the gas 
is isothermal, pgas oc n^^^, then the gas distribution is given by the following 
(Cavaliere & Fusco-Femiano 1976) 

Pgas = Pgasim + {R/ Rc)^]^'"''^ , (8) 

where Re is the core radius and (3 is given by, 



/3 = ^^?^, (9) 



and (T is the line of sight velocity dispersion. Both quantities are derived from 
the observed surface brightness profile which is found to be well characterized 
by a simple analytical form: 

S{x) = So[l + (x/i?,)2)]-3'3,,,+l/2 (^0) 

This functional form gives relatively accurate fits to the data (Jones & For- 
man 1984) except in the central regions of clusters where cooling flows occur. 
Typical values of f3fit^2/3 are smaller than the value obtained using (9) from 
the measurements of Tx and a. This discrepancy is the so called (3 — problem 
and has been thoroughly discussed in the litterature. Some solutions have been 
suggested to solve this problem (Bahcall & Lubin 1994, Evrard 1990, Navarro 
et al. 1995) see also Gerbal et al. (1995). Smaller than this typical values are 
obtained by Durret et al. 1995 with a mean value around 0.4. In their work, 
Durret et al. have analyzed a sample of 12 Einstein clusters with an improved 
method (Gerbal et al. 1994) of analysis which derives the density and temper- 
ature profiles of the X-ray gas by comparing a real cluster X-ray image to a 
"synthetic" image for which the counts predicted to be detected by the IPC was 
calculated by taking into account all the characteristics of the detector such as 
the point spread function, the effective area as a function of radius and energy. 
The ellipticity of the cluster is also taken into account. The resulting simulated 



images are fitted pixel per pixel to observed ones by minimizing the following 
function : 
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A consequence of such flat (small /3) gas density profiles is the derived gas mass 
to dynamical mass ratios (baryon fraction) which are exceedingly large. Another 
interesting result of this analysis is the highly centrally peaked dark matter dis- 
tribution in good agreement with the results based on the imaging and modelling 
of gravitational arcs in clusters (Tyson et al. 1990, Hammer 1991, Mellier et al. 
1993, Wu & Hammer 1993). 
Using isothermality and (8), equation (7) becomes : 



M{r) = 3f3/G 



finip 1 + (r/rcY 



(12) 



with fi = 0.59. This method has been extensively used to derive cluster masses, 
but still one may ask how secure this method is? 

The accuracy of the hydrostatic, isothermal "beta-model" method has been ex- 
amined through hydrodynamical numerical simulations (Schindler et al. 1995, 
Evrard 1996). In particular, Evrard has shown that this method gives remark- 
ably accurate masses inside a radius between 0.5-2.5 Mpc but with a large scatter 
(15 - 30%) (Figure 5). 
However, Bartelmann and Steinmetz (1996) have reached the opposite conclu- 
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Figure 5. Histograms of the estimated mass from the (3 model from Evrard 
1996 

sion, they have found from their gas-dynamical simulations that the (3 — model 

yields systematically low cluster mass estimates . 

Furtheremore, Balland & Blanchard (1995) have discussed the validity of using 
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equation (7) to infer the mass M(r) from the observed temperature T(r). They 
argue that the hydrostatic equihbrium equation is unstable and, using a Monte- 
Carlo procedure, that the resulting accuracy of the mass estimates is rather poor; 
larger than generally claimed. Applying their procedure to the Coma cluster, 
they find a factor of at least 2 uncertainty in the mass inside the Abell radius, 
even when the measurement of the temperature is improved using ROSAT data 
(Figure 6). An alternative way to go round the f3 — model i.e the surface bright- 
ness fitting is not required, has been suggested recently by Evrard (1996). This 
new method exploits an interesting result of his simulations, that is the tight 
relation between the mass and the temperature and uses the resulting scaling 

relations: r^Qo{Tx) oc T^ and M^oq{Tx) oc T^ which lead to more accurate 
masses and the scatter found in the /3 — model is then eliminated (Figure 7). Of 
course such conclusions are given in the frame of numerical simulations which 
simulate clusters in "somehow" perfect conditions. For example their analysis 
uses the clusters emission- weighted temperature which comes from their simu- 
lations and not from cluster spectra. Furthermore, the simulated X-rays images 
can be analyzed out to large radii which is not generally the case in real ob- 
served X-rays ones. Finally, the /3 — model method is based on the assumption 
of spherical symmetry. However, more often clusters exhibit a more complex 
morphologies due to the presence of substructures. Numerical simulations have 
demonstrated that masses of clusters which are undergoing a merging event, are 
generally under-estimated because part of the energy of the gas is in the kinetic 
form due to the bulk motion rather than in the thermal form, therefore the 
temperature of the gas is underestimated and so are the clusters masses. The 
underestimation of the mass due to the presence of substructure can reach 40% 
(Schindler 1996). 
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Figure 6. This plot shows how the mass profile is poorly constrained by 
observed temperature profiles even with the recent results on Coma cluster 
from ROSAT (Briel et al. 1992) 



4.3. The Gravitational lensing mass estimates 

The discovery of giant blue luminous arcs in clusters A370 and CI 2244-02 (Sou- 
cail et al. 1987, Lynds & Petrosian 1989) has provided the first observational 
evidence that clusters of galaxies may act as gravitational lenses on background 
galaxies, a possibility which was first discussed by Noonan (1971). Gravita- 
tional lensing provides a very powerful tool to directly measure the projected 
mass distribution. This method, presents many advantages over the X-ray mass 
estimates, for example, it does not require any assumption on the mass distri- 
bution or on the dynamical state of the cluster. Since the pionnering work by 
Tyson et al. (1990). It has become more and more common to use weak gravi- 
tational lensing to map the dark matter distribution in clusters. Detailed study 
of image formation through gravitational lensing can be found in the review 
by Schneider, Ehlers and Falco ( 1992) and Fort & Mellier (1994). I will just 
summary very briefly the manifestations of the lensing effect and the way the 
lensing masses are derived. The lensing effects can be divided into two main 
regimes depending on the lens configuration : 

1- The strong lensing regime : 

The distorsion of distant galaxies by foreground clusters of galaxies gives rise to 
the spectacular strong arcs observed in the central regions of clusters e.g. A370 
corresponding to a large magnification and strong distorsion. The arclet regime 
is intermediate between the arc and the weak distorsion regimes. 

2- The weak lensing regime or weak shear: 
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mean = 1.00 
s.d. = 0.15 
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Figure 7. Histograms of the estimated mass from the scaUng law from 
Evrard 1996 



The first observational detection using optical galaxies as sources is due to Tyson, 
Valdes and Wenk (1990). In this case, each source produces only one image which 
experiences only a weak distorion of its shape. 

The strong lensing regime constrains the total mass enclosed within the " Ein- 
stein radius" , while weak shear effects determine the distribution of the mass at 
the outer regions (see Brainherd in these proceedings). 

Constraints from Strong Lensing The projected cluster mass within the Ein- 
stein radius rE of arc or arclet can be easily derived if one assumes a spher- 
ical matter distribution for the lensing cluster and assuming that the system 
observer-lens-source is aligned along the line of sight 



M,e 
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(13) 



where Ejens 



Ds 



is the critical mass density with Ds, Di and Dig being 



AttG D,D,s 

the distance to the source (the galaxy), distance between the source and the 
lensing cluster and the distance between the source and the cluster respectively. 
For more complex configurations, cluster masses are estimated by lens modelling 
(see Fort &: Mellier 1994 for a review). This method, however gives the mass 
inside the radius where the arcs are observed which are usually very small ~ 50 
kpc. 



Constraints from weak lensing To construct the surface mass density profile 
one uses the statistic suggested by (Fahlman et al. 1994) 



X{ri,r2) = ^{n) - ^{n <r <r2) 
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where < e > is the mean tangential component of the image ellipticities. This 
method has been successfully applied to several clusters (see Table 1). How 
reliable is the lensing method? The main shortcomings of the lensing method is 
that the application of equation (14) requires an estimation of X^crit ^^^ therefore 
the knowledge of the redshift of the sources which is difficult to obtain. This 
may introduce large uncertainty in the mass especially for distant clusters. On 
the other hand it is not possible to obtain a true value of the mass only from the 
shear map, even in the best case where the sources redshift is known because of 
the degeneracy due to the fact that the addition of a constant mass plane does 
not induce any shear on background galaxies. This degeneracy may be broken 
by measuring the magnification fi of the background which gives an absolute 
measurement of the mass. Broadhurst et al. (1995) have proposed a very nice 
method to measure ^ by comparing the number count in a lensed and unlensed 
field. They find that depending on the slope of the number count in the reference 
field s=dlogN(m)/dm, they observe more or fewer objects in the lensed field. In 
the case where blue galaxies are selected, the counts are unaltered, since the slope 
is in this case equal to the critical value s=0.4. This method has been applied 
successfully to the cluster A1689 by Broadhurst (1995). The weakness of this 
method is that it requires the measurement of the shape, size and magnitude of 
very faint objects. Van Waerbeke et al. (1996) have recently suggested a new 
method to analyzje the lensing effects which avoids the measurement of the shape 
parameter. But still, the weak lensing method leads to very encouraging results 
and promises to yield unambigeous information about the mass distribution in 
the near future. 

4.4. Comparison between X-ray and lensing cluster mass estimates 

Miralda-Escude & Babul (1995) have raised an interesting puzzle. They found 
from their analysis of Abell clusters A2218, A1689 that the mass in the central 
part of the cluster inferred from the strong lensing method is greater than that 
derived from the X - ray method by a factor of 2 - 2.5. Wu & Fang have gathered 
all the clusters for which the mass has been estimated and compared the X-ray 
to lensing masses. They have found a systematic discrepancy between the two 
masses at small radius ~ 0.25 h^Q Mpc which vanishes at larger radii. However 
the lensing and the X-ray information in their sample do not come from the same 
cluster. Early studies based on both optical and lensing observations have led to 
the same conclusion: there is a cluster mass discrepancy by the same factor (Wu 
Fang 1994, Fahlman et al. 1994). But, it seems from a recent statistical analysis 
that viiial masses are consistent with gravitational lensing masses (Wu & Fang 
1997). The disagreement between the lensing masses and X-ray masses may be 
due to the fact that X-ray analysis, namely the (3 — model, underestimates the 
masses. Indeed, the assumption of hydrostatic equation may be invalid, because 
of several reasons, non-thermal pressure, merging effects, a multi-phase medium, 
unstability of the equilibrium equation etc... Unfortunately, it is hard to quantify 
all these effects and to know which is the most important one. 
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Table 1. Lensing masses for a sample of elusters. Ref. (l)Kncib et al. 
1993; (2,3) Tyson & Fisher 1995; (4,5) Miralda-E & Babul 1995, Squires et 
al. 1996a; (6) Kneib et al. 1996,(7) Small et al. 1995b; (8,9) Pello et al. 1991, 
Squires et al. 1996b;(10) Giraud 1988; (11, 12) Wallington & Kochanek 1995, 
Bonnet et al. 1994; (13)Mathez et al. 1992,(14) Hammer et al. 1989, (15) 
Fahlman et al. 1994; (16) Luppino & Kaiser 1996;(17)Smail et al. 1995a; (18) 
Allen et al. 1996; (19)Scliindler et al. 1995 



5. Discussion and conclusion 

Dynamical analysis of clusters of galaxies have led to two important results: 
the presence of large amount of dark matter and the evidence of high baryonic 
fraction, both have implications on cosmology through 0,q and fi^, the density 
of the Universe and its baryon content respectively. Estimating the masses of 
clusters of galaxies, is not straightforward, because it depends on the validity of 
the assumptions underlying the method from which the mass is determined, the 
mass-foUows-light in the case of the virial masses, the hydrostatic equilibrium 
and isothermality of the gas for the X-ray mass determination. Gravitational 
lensing methods provide with a new strong tool to constrain both the amount of 
mass and its distribution. Comparing the X-rays to lensing masses give rise, at 
least in the inner part of the cluster, to the mass discrepancy problem. The most 
probable explanation, would be the underestimation of the X-ray mass. The in- 
teresting implication, is that clusters would be more massive than we think, and 
the ratio of gas mass to total mass (the fraction of baryons) could be in more 
better agreement with nucleosynthesis predictions and an rio=l Universe. Fi- 
nally, thanks to new recent set of observations, it appears that virial masses are 
in good agreement with the lensing masses (Wu & Fang 1997), if this result is 
true, that means that the virial masses are accurate and one may conclude that 
indeed, the mass follows light, since it is only in this case that the virial method 
gives accurate mass determination (Sadat 1995). 
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